Faraday resonance in dynamical bar instability of differentially rotating stars 
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We investigate the nonlinear behaviour of the dynamically unstable rotating star for the bar mode 
by three-dimensional hydrodynamics in Newtonian gravity. We find that an oscillation along the 
rotation axis is induced throughout the growth of the unstable bar mode, and that its characteristic 
frequency is twice as that of the bar mode, which oscillates mainly along the equatorial plane. A 
possibility to observe Faraday resonance in gravitational waves is demonstrated and discussed. 

PACS numbers: 97.10.Kc, 04.25.D-, 04.30.Db, 95.30.Lz 



I. INTRODUCTION 

Parametric resonance is widely observed in hydrody- 
namics, nonlinear optics, chemical reaction and classi- 
cal oscillatory systems. It is also interesting from the 
bifurcation theory and the pattern formation. The sci- 
entific study in the fluid mechanics dates from the ex- 
periments by Faraday in 1831, and is therefore named 
Faraday resonance (e.g., [![). Nonlinear dynamics ex- 
hibits mode interaction of oscillation in different direc- 
tion, and possibly causes the resonant growth of a par- 
ticular mode. Recently, experimental studies of Faraday 
resonance demonstrate that the system of fluid mechan- 
ics [IJ and that of Bose-Einstein condensate [H work per- 
fectly. These agreements in different fields also suggest 
that Faraday resonance may also occur in an astrophys- 
ical context. Quasi-periodic oscillation in gravitational 
waves from dynamical/secular instabilities is expected to 
be excited throughout rotating core collapse, and may 
drive resonant growth. 

Dynamical bar instability in a rotating equilibrium star 
takes place when the ratio (3 (= T/W) between rotational 
kinetic energy T and the gravitational binding energy W 
exceeds the critical value /3dyn (~ 0.27 for an uniformly 
rotating incompressible body in Newtonian gravity d). 
Determining the onset of the dynamical bar-mode insta- 
bility, as well as the subsequent evolution of an unsta- 
ble star, requires a fully nonlinear hydrodynamic simula- 
tion. Simulations performed in Newtonian gravity (e.g., 
I, S 0, 1, I, M M Gl 0, El) have shown that /3 dyn 
depends very weakly on the stiffness of the equation of 
state. /3d yn becomes small for stars with high degree of 
differential rotation lK ll| [ijj • Simulations in relativis- 
tic gravitation [lj|, 

mm 

have shown that /3dyn decreases 
when increases the compactness of the star, indicating 
that relativistic gravitation enhances the bar-mode in- 
stability. Recent numerical simulations show that dy- 
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namical bar instability can occur at si gnifi cantly lower (3 
than the threshold (3 CIt « 0.27 [HI HI M H| • These 
recent findings can be classified into the category of a low 
T/W dynamical instability. This instability may be trig- 
gered by the corotation resonance [III H3| , which is com- 
pletely different from the standard dynamical bar-mode 
instab ility triggered by a certain magnitude of rotation 

Bam. 

Our main concern in this paper is not to determine 
the onset of the instability, but to study the dynamical 
features of the bar. For this purpose, we numerically 
study the growing behaviour of the azimuthal modes in 
the nonlinear regime for a longer timescalc. One inter- 
esting issue of nonlinear evolution is the possibility of 
resonant growth of other azimuthal modes triggered by 
the dynamical bar-mode instability. One candidate for 
such resonance is Faraday resonance, which is excited by 
the external periodic force. According to the linear ap- 
proximation of the velocity potential by using an incom- 
pressible inviscid liquid in a rectangular tank, the time- 
dependent behaviour of the liquid surface is expressed 
by the Mathieu's equation (e.g., [30]). The dynamically 
unstable bar mode may work for other azimuthal oscilla- 
tion modes as an external periodic force. The oscillation 
is not exactly periodic, but rather quasi-periodic, and 
may trigger a parametric resonance. 

The other interesting issue of nonlinear evolution is 
the duration of the bar shape, when it forms. This is 
quite important for gravitational wave detection. We ba- 
sically believe that once the dynamical bar instability 
takes place, the system generates quasi-periodic gravi- 
tational waves for a period sufficient enough to be de- 
tected in the ground-based gravitational wave detectors. 
The only causes to destruct a bar are dissipative effects 
such as viscosity and gravitational radiation. The typ- 
ical timescale of such effects takes place in the secular 
timescale, which is much longer than the dynamical one 
of the system. Therefore the standard picture is that the 
bar can persist in its shape until the secular timescale. 
However, recent numerical simulation shows that a bar 
destructs its shape in the dynamical timescale 20]. The 
authors argue a possible cause of the destruction of bar as 
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azimuthal mode coupling. Although there was in the past 
a debate between the two numerical simulations about 
the persistence of a bar [12, HH , the different outcomes 
were considered as the different accuracy level of the cen- 
ter of mass condition at that time. Since there is the only 
one group that claims the destruction of bar structure in 
the dynamical timescale with a satisfaction of the center 
of mass condition, it is worth investigating the destruc- 
tion of a bar employing a different computational code. 
In order to focus on this topic, it is sufficient to inves- 
tigate this topic in three-dimensional hydrodynamics in 
Newtonian gravity. 

This paper is organized as follows. In Sec.[Tl]we present 
the basic equations of our hydrodynamic simulation in 
Newtonian gravity. In Sec. IIHI are discussed the nu- 
merical results of our findings of Faraday resonance. In 
Sec. IIVI we briefly summarize our findings. Throughout 
this paper, we use the geometrized units with G = c = 1 
[36| and adopt Cartesian coordinates (x, y, z) with the 
coordinate time t. Note that Latin index takes (x, y, z). 



where Sv = 25xd{V l , 5x(— Ax — Ay — Az) is the local 
grid spacing and where we choose the dimensionless pa- 
rameter C v i s = 2. When evolving the above equations we 
limit the stepsize At by an appropriately chosen Courant 
condition. 

The gravitational potential is determined by the Pois- 
son equation 



A$ = 47rp, 
with the outer boundary condition 

M diX 



Here M is the total mass 



+ 0(r- J ). 



M= I pdx A 



(2.7) 



(2.8) 



(2.9) 



and di is the dipole moment 



II. BASIC EQUATIONS 

A. Newtonian Hydrodynamics 

We construct a three dimensional Newtonian hydro- 
dynamics code assuming an adiabatic T-law equation of 
state 



P = (T - l)pe, 



(2.1) 



where P is the pressure, T the adiabatic index, p the mass 
density and e the specific internal energy density. For 
perfect fluids the Newtonian equations of hydrodynamics 
consist of the continuity equation 



dp d(pv l ) 
dt dx l 



0. 



the energy equation 

de d(evi) 
dt dxi 

and the Euler equation 

d(pvi) d(pViV j ) 



.I e -(r-i)p. ^1 

r vls dx* 



(2.2) 



(2.3) 



dt 



dxi 



d{P + P vis ) 

dx % 



<9$ 

"5?- (2 - 4) 



Here v % is the fluid velocity and $ the gravitational po- 
tential; and e is defined according to 

e={pef' v . (2.5) 

We compute the artificial viscosity pressure P v i s from [32| 

_ / C vis p(Svf, for^<0; 



0, 



for Sv > 0, 



(2.6) 



di = / pxidx 
'v 



B. Initial Data 



(2.10) 



As initial data, we construct differentially rotating 
equilibrium models with an algorithm based on Hachisu 
[331 ] and adopt cylindrical two dimensional coordinate to 
compute the axisymmetric equilibrium star. Individual 
models are parameterized by the ratio of the polar to 
equatorial radius R p / R oq , and a parameter of dimension 
length d that determines the degree of differential rota- 
tion through 



n 



Jo 



d 2 



(2.11) 



Here Q is the angular velocity, jo a constant param- 
eter with units of specific angular momentum, and w 
the cylindrical radius. The parameter d determines the 
length scale over which Q changes; uniform rotation is 
achieved in the limit d — ► oo. For the construction of 
initial data we also assume a polytropic equation of state 



P 



Kp 



1+1/n 



(2.12) 



where n = 1/(T — 1) is the polytropic index and k a 
constant. In the absence of shocks, the polytropic form of 
the equation of state is conserved by the T-law equation 
of state (Eq. [2~T] ). 

We also compute the virial identity, which is identically 
zero in the equilibrium star, to show the accuracy level 
as 



\2T tl 



-W + 3U\ 
~W 



(2.13) 
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where 

T to t = \l pv l Vl d 3 x, (2.14) 

w = -\ J P®d 3 x, (2.15) 

n = J Pd 3 x. (2.16) 

Note that we have divided the value by a gravitational 
binding energy W so that the value Vn w t is regarded as 
a relative error of the system. We summarize our four 
different rotating equilibrium stars in Table |TJ 



in equations (|2~19l) - (|2~22"1) (see Ref. [H|). 

The spectrum of gravitational waveform can be com- 
puted as 

S (x,z) = \~ h ^\ 2 + \h^' z) \ 2 , (2.24) 

where 

= f dth^$e iut . (2.25) 



D. Diagnostics 



C. Gravitational Waveforms 

We compute approximate gravitational waveforms by 
evaluating the quadrupole formula. In the radiation zone, 
gravitational waves can be described by a transverse- 
traceless, perturbed metric hj^ with respect to a flat 
spacetime. In the quadrupole formula, hf T is found from 



■TT 



2d?_ 
rdt 2 ^ 



TT 



(2.17) 



where r is the distance to the source, Jy the quadrupole 
moment of the mass distribution (see Eq. [36.42b] in 
Ref. [3ij]), and where TT denotes the transverse-traceless 
projection. Choosing the direction of the wave propaga- 
tion to be along the z-axis (rotational axis of the equi- 
librium star), we determine the two polarization modes 
of gravitational waves from 



h^ = \{h T ^-h T ^) and h!£> = h™. 



For observers along the z-axis, we thus have 



rh, 



M 



rh 



M 



1 d 
2Mdt 
1 d ■ 

Mdi 



Ixy 



(2.18) 

(2.19) 
(2.20) 



Using the same procedure, the observers along the x-axis 
detect the wave propagates as 



,(1) 
rh + 

M 



1 d , • 
2~M~dt 



(lyy Izz), 



rh 



(x) 



M 



- Mdt vz - 



(2.21) 
(2.22) 



The number of time derivatives that have to be taken 
out can be reduced by using the continuity equation 
(Eq. E3) 



Iij — 



{pv l x i + px l v j )d 3 x, 



(2.23) 



We monitor the conservation of mass M (Eq. |2.9| ). 
angular momentum J 



J 



p(xv y - yv x )d x, 



and the location of the center of mass x 



CM 



CM 



px L d 3 x. 



(2.26) 



(2.27) 



Due to our flux-conserving difference scheme the mass M 
is also conserved up to a round-off error, except if matter 
leaves the computational grid. 

To monitor the development of the azimuthal modes 
(to = 1, 2, 3, 4) and the one in the z-direction, we com- 
pute the following five diagnostics 



D = (e lmip ) , 

\ / 771=1 

If x + iy 3 
M J ^fx 2 + y 2 

1 f {x 2 - y 2 ) + i(2xy) 3 



M J x 2 + y 2 



d 3 x, 



(2.28) 



(2.29) 



O = (e tmip ) 

\ I 7 



1 f x(x 2 - iy 2 ) + iy(3x 2 - y 2 ) 
P — 



M J r (x 2 + y 2 f/ 2 

M 4 = (e lmip ) . 



d 3 x, (2.30) 



1 f (x 4 -6x 2 y 2 + y 4 )+i(4x 2 y 2 (x 2 - y 2 )) 



M 



9- 



(x 2 + y 2 ) 2 



d 3 x, 



D, 



MR r , 



p\z\d 3 x, 



(2.31) 
(2.32) 



where a bracket denotes the density weighted average. 
When we compute the four diagnostics in the equatorial 
plane (D^, Q (cq \ (oq) , Afi° q) ), we change the integral 
volume from d 3 x to dxdy and M to M eq . Note that M oq 
is the rest-mass density integrated only in the equatorial 
plane. 
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FIG. 1: Diagnostics \D\, \Q\, \0\, \M^\ as a function of t/P c for four different rotating stars (see Table [I]). Solid, dashed, 
dotted, and dash-dotted lines denote \D\, \Q\, \0\ and \Mn\, respectively. We terminate our simulation when the relative error 
of the rest mass exceeds ~ 0.01% for models I and II. Hereafter P c represents the central rotation period at t — 0. 



We also compute the spectra of the above five diagnos- 
tics as 



1 

M 



dt / d xpe 



MIL 



dt / d 3 x\z\(p - p avg )e 



(m= l,---,4), 
(2.33) 
(2.34) 



Since D z does not oscillate around zero as we show later, 
we have subtracted the time averaged density p avg from 
the original one to compute the spectrum. Note also that 
we have only integrated D z in time after its first global 
maximum. 



J 



P = Poq 



2xy - y 2 



III. NUMERICAL RESULTS 

Here we show our evolution of the differentially rotat- 
ing stars. We terminate the integration when the relative 
error of the rest mass exceeds ~ 10 , since the only vi- 
olation of the rest mass conservation is caused by the 
matter outflow at the outer boundary of computation. 
We also terminate the integration when the time exceeds 
20 ~ 40 central rotation periods, which are sufficient to 
enhance all m modes. Note that our code never crashes 
throughout the evolution. 

To enhance any dynamically unstable mode, we dis- 
turb the initial equilibrium density p ea by a non- 
axisymmetric perturbation according to [37[ 



R 2 



aw 



x 4 — 6x 2 y 2 + y 4 + Axy(x 2 — y 2 ) 



R 4 



(3.1) 



where we set 8^ = 8^ = 10 2 . 

We study four different differentially rotating stars, 
which are detailed in Table [J to investigate the nonlinear 
behaviour of the non-axisymmetric dynamical bar insta- 
bilities. We choose the axis of rotation to align with the 
z axis, and assume planar symmetry across the equator. 
We choose the Cartesian coordinates with the computa- 
tional grid points 401 x 401 x 101 covering the equatorial 



diameter of the equilibrium star as 121 grid points. 

We show the amplitudes of our four diagnostics for all 
four models in Fig. [TJ At the first stage of evolution, 
the m — 2 diagnostic grows exponentially in models I, 
II, and III, while it stays around the amplitude of t = 
in model IV. Therefore the star of models I, II, and III 
is determined as dynamically unstable against bar mode, 
while that of model IV is stable. For the dynamically bar 
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FIG. 2: Center of mass as a function of t/P c . Solid, dashed line denotes the one of x and y direction, respectively. 
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FIG. 3: Linear momentum as a function of t/P c . Solid and dashed line denotes x and y component of the linear momentum, 
respectively. M and fi e represent the total rest-mass and the angular velocity at the equatorial surface at t — 0, respectively. 



unstable stars (models I, II, and III), the m = 2 diagnos- 
tic grows exponentially but the other remaining m modes 
do not grow at the first evolution stage when imposing a 
small perturbation (Fig. QJ. This result is consistent to 
the linear perturbation analysis of the dynamically bar 
unstable stars, which shows that the only dynamically 
unstable m mode is m = 2. Also the result confirms us 
that the amplitude of perturbation at t = is adequate to 
treat the system linearly (S 1 - 2 ^ w 10~ 2 ). After that stage 



the m = 4 diagnostic grows exponentially because of the 
secondary harmonic of m — 2 mode, and then the odd m 
modes are also enhanced. The odd m modes are excited 
even if we do not impose the perturbation of their corre- 
sponding modes in the equilibrium star, since the finite 
differencing scheme always generates a small amount of 
all m modes (Fig. [T]) . However a small fluctuation at the 
wavefront should occur in nature so that the existence of 
all m modes, when the bar forms, are quite natural in 



6 




TABLE I: Four different rotating equilibrium stars in New- 
tonian gravity of F — 2, d/R cq = 1. 

Model R p /R cq T/W y Ntw 

I 0.225 0.281 8.29 x 10 -5 

II 0.250 0.277 8.79 x 10~ 5 

III 0.275 0.268 7.95 x 10~ 5 

IV 0.300 0.256 9.47 x 10~ 5 



reality. 

We have monitored the center of mass and the linear 
momentum throughout the evolution to guarantee that 
we do not impose any additional physics in the system. 
Fig. [5] shows the center of mass of the four different stars 
throughout the evolution. We have confirmed that the 
numerical error only allows the star to change the cen- 
ter of mass within the one computational grid. We have 
also checked the linear momentum conservation in Fig. [31 
which shows that the relative error is less than 1% of the 
total value constructed by the total mass and the ve- 
locity at the equilibrium equatorial surface. In order to 
check whether the center of mass condition significantly 
affects the diagnostics, we have also computed the follow- 
ing two types of diagnostics in the equatorial plane (e.g. 
[251). One is the diagnostic with the same coordinate 
as in the simulation, while the other is the one with the 
coordinate where the center of mass is adjusted to zero 
in every snapshots. Since the equatorial diagnostic (Fig. 
SI) reproduces all characteristics of the one obtained from 
the three-dimensional computation (Fig. [1]) , the equato- 
rial diagnostic may represent the three-dimensional one. 



We compare Figs. @] and to focus on the effect of the 
center of mass condition on the diagnostics. For models 
I and II, the adjustment of the center of mass reduces 
the amplitude of Z)( oq ' and 0^ eq 'for t < 5P C . However 
the condition does not change the exponential growth of 
£>(° q ) and 0(° q ) after t > 5P C . For models III and IV the 
amplitude of £)( cq ) has been reduced so that the system 
is stable to m — 1. Therefore the linear growth of £)( eq ) 
in models III and IV is the outcome of the violation of 
the center of mass condition. 

We also show our equatorial and the meridional den- 
sity snapshots throughout our integration in Figs. [6] and 
[7] The symmetry breaking of the bar structure occurs 
clearly at the time when the spiral arm forms in the 
equatorial snapshot and in the meridional plane. This 
becomes clear when we focus on the final snapshots of 
models I, II, and III. 

We also in Fig. [5] show the diagnostics which contain 
both amplitude and phase. In order to make the picture 
clear, we first concentrate on the model III, the weakest 
dynamically unstable bar system of three models. 

The behaviours in the diagnostics are clearly under- 
stood once we compute the spectra of the diagnostics 
(Fig. [5]) . From the spectra we find the following two re- 
markable issues. One is that the spectra |-Fi| 2 , I-F2I 2 , 
|F 3 | 2 take a peak around Wbar ~ 5 ~ 6P ( T 1 for models 
I, II, III, and the other is that |F 3 | 2 , |F 4 | 2 , \F Z \ 2 take a 
peak around ix> qU ad ~ 2ubar ~ 10 - 12P c _1 for bar unstable 
stars. Combining the present feature with the behaviour 
of the five diagnostics explained before (Fig. [1]) , the dy- 
namically unstable bar acts as follows. 

Firstly the m = 2 mode grows and acts as a domi- 
nant mode of all because of the dynamical bar instabil- 
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FIG. 5: Same as Fig. [4] but the center of 
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is adjusted to zero in every snapshots. 



ity. Next the to = 4 mode grows because of the secondary 
harmonic of the m = 2 mode. In fact the saturation am- 
plitude of the to = 4 is approximately »s 0.3 for model 
I, 0.2 for model II, and 0.04 for model III, all of which 
are the order of the square of the saturation amplitude 
of the m = 2 (~ 0.6 2 for model I, 0.5 2 for model II, 
0.2 2 for model III). After that Faraday resonance occurs, 
which is clearly found in both D z and P 2 | 2 from the fact 

^quad ~ 2wbar- 

Note that Faraday resonance occurs in the fluid me- 
chanics when the oscillation of the vertical direction is 
twice (2w) as much as the one in the horizontal direction 
(uj) in the weakly nonlinear interaction [3,[3(|. The rea- 
son why the resonance does not clearly appear in model 
I is either the strongly nonlinear effect or the insufficient 
duration time of quasi-periodic oscillation for computing 
the spectrum. Then, there is a resonance between m = 1 
and to = 2, to = 3 and to = 4. The possibility of such 
resonances is three wave interaction: either to = 1 (wbar) 
and to = 2 (t^bar) generates to = 3 (wbar+^bar) or to = 3 
(2wbar) and to = 2 (uw) generates to = 1 (2u> baT — u> baI ) 
in the dominant part. It is the fact found in the nonlinear 
behaviour of the dynamically unstable bar system. 

The gravitational waveform and its spectrum have 
been computed by the quadrupole formula observed 
along the rotational axis and in the equatorial plane 
(Figs. [10] - [T3|) . There are two remarkable features in 
gravitational waves from the viewpoint of nonlinear be- 
haviour. One is that the quasi-periodic oscillation does 
not last until the radiation reaction timescale but de- 
cays because of the symmetry breaking of the dynamical 
bar. The duration period is related to the degree of non- 
linearity of the bar mode instability, which is estimated 
from the inclination angle of the amplitude of the m = 1 



(Rc[D]) and to — 3 (Re[0]) diagnostics. In the present 
case, the duration period of the bar structure is estimated 
as ~ 10P C for model I, ~ 15P C for model II, and ~ 35P C 
for model III. The other is that Faraday resonance has 
clearly appeared in the spectrum of gravitational wave- 
form observed at least in the equatorial plane. Since 
we adopt quadrupole formula to compute gravitational 
waves, the higher order harmonics of the unstable bar 
mode such as m = 4 mode cannot be seen in this spec- 
trum. Therefore a peak around u> w 12P ( 7 1 in Fig. [TTJ 
indicates the fact of an oscillation along the z-axis, which 
is the evidence of Faraday resonance. We have also com- 
puted the gravitational waveform and its spectrum ob- 
served along the z-axis and found that there is no peak 
around u « 12P C -1 in model III (Fig.UJ). The fact also 
supports that a peak around id k 12P C _1 in Fig.[TTJis the 
outcome of Faraday resonance, since an oscillation along 
z-direction can be clearly observed by gravitational waves 
in the equatorial plane, not in the rotation axis. When 
we increase the degree of nonlinearity, the above feature 
of the Faraday resonance in gravitational waves can be 
also seen in the equatorial plane. 

IV. CONCLUSION 

We investigate the nonlinear effects of dynamically bar 
unstable stars by means of three dimensional hydrody- 
namic simulations in Newtonian gravity. In order to fol- 
low the bar shape as long as possible, the initial ampli- 
tudes for odd azimuthal perturbations are significantly 
suppressed in our models. 

We find interesting mode coupling in the dynamically 
unstable system in the nonlinear regime, and that only 
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FIG. 6: Density contours in the equatorial plane throughout the evolution. Snapshots are plotted at (t/P c , pmax/pmax, d) = 1(a) 
(3.70xl(T 4 , 1.00, 0.25000), 1(b) (2.40, 1.21, 0.25000), 1(c) (5.92, 1.27, 0.25000), 1(d) (8.87, 1.46, 0.25000), 11(a) (3.29xl0~ 4 , 1.00, 
0.20625), 11(b) (3.95, 1.20, 0.20625), 11(c) (7.90, 1.29, 0.20625), 11(d) (11.85, 1.58, 0.20625), 111(a) (2.98 x 10~ 4 , 1.00, 0.25000), 
111(b) (11.93, 1.13, 0.25000), III(c) (23.85, 1.19, 0.25000), 111(d) (35.78, 1.27, 0.25000), IV(a) (2.75 x 10~ 4 , 1.00, 0.25000), IV(b) 
(12.10, 1.05, 0.25000), IV(c) (24.20, 1.10, 0.25000), IV(d) (36.30, 1.16, 0.25000), where p max is the maximum rest mass density 
and Pmix is the maximum rest mass density at t = 0. The contour line denotes p/p max = 10 ( - 16 ~ I - )d (i = 1, . . . , 15). Hereafter R 
denotes the equatorial radius at t = 0. 



before the destruction of the bar. The quasi-periodic os- 
cillation mainly along the rotational axis is induced. The 
characteristic frequency is twice as big as that of the dy- 
namically unstable bar mode. This feature is quite anal- 
ogous to the Faraday resonance. Although our finding is 
only supported by the weakly nonlinear theory of fluid 
mechanics, we have also found the same feature of para- 
metric resonance even in the strongly nonlinear regime. 
There is one qualitative difference between Faraday res- 



onance and our numerical result. Faraday resonance has 
lower frequency than that of the forced oscillation, while 
our result has higher frequency than that of the bar un- 
stable mode. The fact can be understood by the differ- 
ent regime of the media. Since the media of the rotating 
star is a perfect fluid, which is only contained inside the 
star, there should be a cutoff frequency to be amplified. 
In fact, introducing a cutoff frequency with a polar ra- 
dius of the star and the sound speed computed by the 
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FIG. 7: Same as Fig. [5] but in the meridional plane. 
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FIG. 8: Diagnostics 5ft[-D], 5R[0], 5ft[M4], and D z as a function of t/P c for four different rotating stars (see Table|I]). 



mean density, the cutoff frequency w cu t is estimated as 
Lu cut ~ Wbar- Therefore a higher frequency than the bar 
is amplified. The fact also indicates that our finding is 
interpreted as a parametric resonance. 



We also discuss the detectability of Faraday resonance 
by gravitational waves. The characteristic frequency and 
the amplitude of gravitational waves can be estimated as 
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FIG. 9: Spectra \F m \ 2 and \F Z \ 2 as a function of ujP c for four different rotating stars (see Table|I]). Solid, dashed, dotted, and 
dash-dotted line of \Fm\ 2 denote the values of m = 1, 2, 3, and 4, respectively. 




FIG. 10: Gravitational waveform (+ mode) observed at the x axes for four different rotating stars (See Table Note that x 
mode is identically zero because we adopt the equatorial symmetry. 
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FIG. 11: Spectra of gravitational waveform (+ mode) observed at the x axes for four different rotating stars (See Table|l]). 
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FIG. 12: Gravitational waveform observed at the z axes for four different rotating stars (See Table H}. Solid and dashed line 
denotes + mode and x mode, respectively. 



where d b s is the distance from the observer. If the bar 
formation occurs in Virgo cluster, quasi-periodic wave- 
form can be detected in the second generation of grav- 
itational wave detectors such as Advanced LIGO, Ad- 



vanced Virgo, and Large-scale Cryogenic Gravitational 
wave Telescope (LCGT), or in the third generation Eu- 
ropean Gravitational Wave Observatory. The frequency 
of the parametric resonance is around twice as big of that 
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of the bar unstable mode, and the amplitude of the para- 
metric resonance is roughly two orders lower (ps 1%) then 
that of bar unstable mode. The detection of gravitational 
waves from parametric resonance may explore the non- 
linear phase of the dynamically bar unstable stars such 
as determining the saturation amplitude of gravitational 
waveform from the bar unstable system, parametric res- 
onance, and the duration period of the bar structure. 
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